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1.  INTRODUCTION 


The  purpose  of  this  report  is  to  provide  a method  for  calculating  the 


scattered  electromagnetic  field  produced  by  illumination  of  a rough  surface.  At 


of  illumination 


polarization.  Over  a substantial  portion  of  the  report,  mostly  for  purposes  of 


comparison  with  other  work,  the  incident  field  will  be  taken  to  be  a plane  wave 


with  simple  polarization  (usually  IM) . When  periodicity  is  assumed  for  the 


surface,  the  incidence  angle  of  the  illuminating  radiation  will  be  restricted  to 


certain  angles  in  order  to  achieve  an  overall  periodicity  compatible  with  the 


when  given  explicitly  by  z = f(x,y)  we  allow  -“>  * x * 00 


vacuum  (z  > f(x,y),  index  of  refraction  equal  to  one)  and  a perfect  conductor 


(z  < f(x,y)  conductivity  infinite).  Finally,  we  assume  that  f(x,y)  and  its  first 


and  second  derivatives  are  continuous.  Where  necessary  for  comparison  with  other 


work 


When  periodicity 


with  the  attendant  requirement  that 


in  f (x,y)  is  assumed  the  Fourier  integral  becomes  the  usual  Fourier  sum.  Note 


that  no  assumption  is  made  concerning  the  smallness  of  the  surface  variation  on 
slopes.  The  only  surfaces  excluded  are  those  having  sharp  points  or  ridges  where 


Formulae  are  given  for  the  calculation  of  the  scattered  field 


Indeed,  without  this  it  is  not  possible  to 


discuss  the  full  spectrum  of  the  scattered  field  in  the  "valleys"  of  the  surface. 


i.e.,  in  the  region  where  z lies  between  maximum  and  minimum  f(x,y),  (0 


f . < z < f ).  Upon  taking  the  asymptotic  form  of  certain  Green's  functions 

min  max 

appearing  in  the  expression  for  the  scattered  field,  the  somewhat  simpler  forms 
of  the  scattered  field  remote  from  the  surface  are  easily  extracted. 

The  general  method  used  in  this  work  is  the  exact  theory  for  the 
reflection  problem  provided  by  a two  stage  calculation.  First,  one  solves  the 
appropriate  integral  equation  system  for  the  surface  currents  induced  by  the 
incident  field.  Second,  these  currents  used  as  sources  furnish  the  structure 
of  the  scattered  field  by  means  of  the  usual  Green's  function  techniques.  In  this 
way  it  is  possible  to  analyze  the  conventional,  or  physical  optics,  approach  to 
the  scattering  problem  and  determine  the  range  of  its  validity.  It  will  be  found 
that  even  apart  from  the  restriction  of  small  surface  roughness,  there  is  a basic 
assumption  of  the  physical  optics  approach  which  is  invalid.  In  order  to  examine 
this  point  further  we  turn  to  a very  brief  recapitulation  of  the  physical  optics 
approach  in  simple  enough  form  to  exhibit  the  difficulties. 


2.  THE  PHYSICAL  OPTICS  SOLUTION 

The  method  presented  in  this  section  is  perhaps  best  exemplified  by 

the  work  of  Rice1  although  the  basic  idea  is  much  older;  the  same  approach, 

, 2 

used  for  acoustical  scattering,  goes  back  at  least  to  Lord  Rayleigh.  Variants 

on  Rice's  original  presentation  have  been  attempted  in  an  effort  to  relax  some 

3 4 

of  the  basic  conditions.  For  example,  Barrick  has  applied  Leontovich 
boundary  conditions  in  order  to  move  away  from  the  restrictions  of  perfect 
conductivity  of  the  surface.  It  may  be  noted  here  that  the  literature  on  this 
problem  is  very  extensive  and  that  no  review  of  the  literature  will  be  under- 
taken in  this  report. 


2. 


In  our  example  we  consider  the  surface  given  by  z = f (x)  with  the 


normal , 


n = n(x)  = (1  + f^r5*  + e2] 


(1) 


f = df/dx;  e , e , and  e are  unit  vectors  in  the  x,  y,  and  z directions 
x x y z 


respectively.  For  the  TM  mode  (E^  = 0)  the  boundary  condition  of  no  tangential 


component  of  electric  field  on  the  surface  assumes  the  simple  form. 


nE  -nE  =0  on  z = f(x) 
Z X X z 


(2) 


which,  upon  retention  of  first  order  terms  assuming  f is  small,  becomes 


E + f E = 0 
X X z 


on  z = f (x) 


(3) 


The  incident  field  has  only  a y component, 


£inc  _ + eik[x  cos  0O  - z sin  0O]  . H £ $inc 

o y o y 


(4) 


It  may  be  noted  that  0 is  measured  from  the  horizontal.  The  corresponding 

o 


electric  field  is  easily  found: 


.me 


H sin  0 e 
o o 


ik  [x  cos  0O  - z sin  0O] 


k o dz 


inc 


,inc  _ 


H cos  0 eik  [x  cos  6°  - z sin  0O]  = i H M 
o o k o 3x 


inc 


(5) 


4 Vot* 

A time  dependence  factor,  e , has  been  assumed  for  all  field  variables 
where  k = 2tt/X  and  X is  the  wavelength  of  the  incident  radiation. 

For  the  scattered  field  the  physical  optics  approach  assumes  a 
superposition  of  upward  propagating  waves,  each  of  which  is  in  form  similar  to 
that  suggested  by  eg.  (4).  The  function  corresponsing  to  4>inc  is 


. s t ik[xp  + z/1  - p^] 

$ (x,z)  = J dp  A(p ) e 


(z  > 0)  (6) 


For  0 < | p | < 1,  these  are  waves  scattered  upward  at  some  fixed  angle  with 
amplitude  proportional  to  A(p).  For  |p|  > 1,  the  scattered  field  consists  of 
trapped  waves  on  the  surface  decaying  in  amplitude  exponentially  with  increasing 


positive  2. 

The  procedure  is  now  straight  forward.  The  scattered  electric  field 
s s 

components,  E^  and  E^,  are  calculated  from  equations  similar  to  eqs.  (5) 

using  $S  in  place  of  4’^nc.  Then  forming  the  total  field 

tot  „inc  j s tot  inc  s 

E - E +E;E  = E +E 
x x x z z z 

the  boundary  conditions,  eq.  (3)  are  applied.  Since  f(x)  as  well  as  its  first 

derivative  is  assumed  small,  an  expansion  in  powers  of  f and  f^  is  made,  and  a 

corresponding  perturbation  expansion  of  A (p ) is  made,  i.e.,  A(p)  = A (p)  + A,  (p)  + 

o 1 

...  Aq(p)  represents  the  contribution  of  a flat  surface  (f(x)  = 0)  while  A^fp)  is 
the  contribution  from  first  order  terms  in  f (x)  and  fx(x)* 


The  zeroth  order  yields  simply. 


A (p)  = 6(p  - cos  6 ) 
o o 


where  6 (x)  is  the  usual  Dirac  delta  function.  The  corresponding  equation  for 


first  order  term  is, 


/ dp  /1_-  p*  A1(p)  elkxp  = - 2[ikf  sin2  0Q  - fx  cos  6q]  elkX  COS  6°  (8) 

- 00 

so  that  an  explicit  solution  for  A^(p)  follows  from  the  Fourier  expansion  of  f(x). 
Suppose  f (x)  has  a horizontal  scale  A and  vertical  scale  a so  that, 

f(x)  = 0/  d<  C(<)  (9) 

.00 

then  with  the  definitions  a = kA,  8 = kA,  and  e ■ o/A  eq.  (7)  yields, 

A (p)  = e c(a(p  - cos  6 ))  [2ia2(p  cos  6 - 1)]  (1  - p2)-*5  d°) 

1 o o 

The  smallness  of  A^(p)  in  comparison  with  Aq(p)  is  indicated  by  the  factor  e in 
eq.  (10).  When  restriction  to  one-dimensional  surface  variation  is  made,  the 
above  result  is  that  of  Rice.  Extension  to  second  order  is  straightforward  and 
will  not  be  repeated  here. 

The  important  point  to  note  is  the  unwarranted  nature  of  the  assumption 

of  upward  propagating  waves  only  (eq.  (6))  coupled  with  application  of  the 

boundary  condition  in  a region  (f  . z $ f ) where  the  whole  spectrum  cannot 

min  max 

be  represented  by  such  a hypothesis.  The  physical  optics  approach  provides  no 
method  by  which  the  requisite  downward  propagating  waves  can  be  introduced.  It 
is  one  of  the  purposes  of  this  report  to  show  explicitly  how  these  downward 
propagating  components  make  their  appearance.  The  successive  order  perturbation 
terms  of  the  physical  optics  approach  can  be  interpreted  as  a form  of  multiple 
scattering  contributions,  but  the  interpretation  is  suggestive  more  than  real 


since  the  successive  terms  correspond  simply  to  iteration  of  a linearized 
(in  e)  integral  equation. 

3.  THE  INTEGRAL  EQUATIONS  FOR  SURFACE  MAGNETIC  FIELDS 

Rather  than  make  any  assumptions  about  the  nature  of  the  scattered 
field,  and  consequently  the  total  field  on  the  surface,  the  viewpoint  taken  here 
emphasizes  the  direct  calculation  of  these  quantities.  The  basic  formula  for 
this  procedure  is  the  Chu-Stratton  integral  equation  for  the  surface  magnetic 
field.  The  second  stage  of  the  calculation  uses  these  surface  fields  as  sources 
for  the  application  of  a suitable  Green's  function  technique  to  yield  the  scattered 
field.  The  integral  equation  for  a perfect  conductor  is 


H(x')  = 2HeXt(x')  + ik3/da  [n(x)  x H(x)J  x ?(x,x' ) eikr 

**  (krr 


(11) 


where  x symbolizes  the  vector  with  components  (x,y,f(x,y))  and  similarly  for 

-►  -V  ->  -►  t 

x',  r(x,  x')  = x - x ' , and  the  integration  is  over  the  whole  surface.  The  vector 

n da  has  components  (-f  , -f  , 1)  dx  dy.  From  the  surface  boundary  condition, 

x y 

■+  -► 

already  used  in  eq.  (11) , n • H = 0 it  is  evident  that  the  third  component  of  the 
magnetic  field  is  determined  in  terms  of  and  so  that  only  two  coupled 

integral  equations  for  and  suffice  to  give  the  whole  surface  field. 

-►ext  • 

H (x)  is  the  incident  magnetic  field  which,  for  plane  wave  illumina- 
tion, has  x and  y components. 


(H®Xt,  Hext)  = HQ(sin  0Q,  -ro)  exp  [ik(x  cos  0o  - f(x,y)  sin  0o] 


(12) 


with  r 


a complex  number  dependent  on  polarization  conditions; 


H is  the  scale 


of  the  incident  field.  Some  possibilities  are:  r ■ 0 gives  TE  mode, 

= ±i  gives  right  (left)  circular  polarization,  while  the  replacement  of 

(sin  0 , -T  ) by  (0,  -1)  gives  the  TM  mode, 
o o 

The  use  of  eq.  (11)  is  much  simplified  by  the  introduction  of  scaled 

variables.  Let  A and  a be  horizontal  and  vertical  scales  respectively  of  f(x,y). 

So  let  x = AC  and  y = An  , and  set  f(x,y)  = ox(£,n).  The  derivatives  of  f(x,y) 

scale  according  to  f (x,y)  = ex  (C,n);  f (x,y)  = ex  (£,n);  f (x,y)  = e2x_  (C,n), 

x s y n xx  s s 

etc. , where  e = 0/a.  The  same  notation  as  follows  eq.  (9)  in  relating  these 
scale  factors  to  incident  wavelength  is  also  used  so  a = kA,  and  0 = kc.  The 
magnetic  fields  are  scaled  by  Hq  so  on  the  surface. 


(H  (x,y,  f (x,y) ) , H (x,y,  f (x,y))  = H (<(>(£, n),  i|>_(£,n)) 
x y o L z 


,„ext , ext,  ..  „ , . . „ . ia[c  cos  0 -ex(£,n)  sin  0 1 

(H  (x,y,  f (x,y) ) , H (x,y,  f(x,y))  - H (sin  0 , -T  ) e o oJ 

x y o o o 

(13) 

The  expression  [(ikr  - 1)  exp  ikr]  /(kr)3  is  replaced  by  (-i^/iOH^^  (aP0)/(aP0) 3^2 
where,  with  x'  = x(C,  h'), 


= [(C-5')2  + (n-n')2  + t2(x-x«)2]3s 


and  h',  is  the  Hankel  function  of  the  first  kind  of  order  3/2. 

j/  £ 

Introducing  the  column  vector  iMC,n)  whose  components  are  ^(£,n)  and 


^2(C»h)  allows  eq.  (11)  to  be  written 


(1) 

• “ H;,'(oiP  ) 


) = 2 »ext(g,,n')  - i — : £g3  M / dn-3/-i.  - ° M(c '.n1  ;S,n)  MS,n>  a5) 

n./nZ  M \3r  z 


2^2*  -<*  -«  (apQ )' 


A similar  column  vector  has  been  introduced  for  the  incident  field 


. ext... 

<P  {£\n’) 


sin  6 


o ia[C  cos  9 - ei(£,n)  sin  9 ] 


,ext  ia[C  cos  0 - et(£,n)  sin  9 ] 

V e o o 

o 


and  M(£'/n'j£/n)  is  a 2 x 2 matrix  with  the  structure, 


MU'.n'i  ?,n)  = Mo(£’,n';  5,n)  + e £,n) 


where 


Mq(C  .n'.C.n)  = 


t - x'  - (n  - n' ) (n  - n' )t^ 
(C  - 5'H  t - x'  - U - ?’)tc 


Tr  TFTn 

M U'.n';  S,n)  = (t  - t')  ^ ^ = (t  - t’)  N2(e,n) 

Z TrT  T 2 

5 n n 


It  may  be  noted  that  e enters  into  eq.  (15)  in  a more  complicated  way  than  the 
usual  Fredholm  equation  of  the  second  kind  since  e also  appears  in  a non-linear 
way  in  pq.  For  this  reason  simple  iteration  of  eq.  (15)  does  not  yield  an 
expansion  of  iM£,n)  in  powers  of  e. 


4.  THE  SCATTERED  FIELD 

Above  the  surface  (z  > f(x,y))  the  total  magnetic  field  is  expressed  in 


non-scaled  coordinates  by 


H(x’)  = H*Xt(x')  - i — i—  k3  / da  [£<x)  x H(x)]  x?(x,x') 


H^i(kr) 


where  x'  now  refers  to  an  arbitrary  point  off  the  surface  (coordinates:  x' , y' , z' ) 


= [(x  - x')2+  (y  - y')2+  (f  - z')2]*5 


Eqs . (18)  and  (19)  are  easily  expressed  in  scaled  coordinates.  From  eq.  (18)  the 
scattered  field  is  readily  identified  as  the  surface  integral  term  with  horizontal 
components  given  in  two  component  column  vector  form  by. 


« Hi^(ap') 


H*(x')  = - i -2-  a / d € / dn  — ms(x';  C,n)  *(£,n) 


4/2tT  -<*  -°°  , ,•>3/2 

(apo> 


in  which  from  eq.  (19), 


, p,_  kx' .2  , ky' .2  . ,kz'  .2-i^ 

po = Cu  - + (n  - a ) + (~  - ET)  ] 


and  the  2x2  matrix  MS(x';£,n)  has  the  structure. 


us . , . _ kz'  _ kz'  1 0 

Mo(x  ' C’n)  = - “ 1 5T  0 1 


M1(x';  £,n)  = 


t - (n  - ^-)t  (n  - 

ah  at, 


kx' . , _ kx' . 

a - T - " - — )T? 


m2(x' ; £,n)  = 


-^N2(C,h) 


(N2(^,n)  is  as  defined  implicitly  in  eq.  (17)) 


For  the  vertical  component  of  the  scattered  field  eq.  (18)  gives, 


1 


H^(x' ) = i — £_  ct3  / d£  J dn  37 2.  ^J~  [v^x' ; C,n)  ^(C.n)  + V2(x' ; C.n)  *2(C.n)] 


4/2tT  — OO  — 00  (ap  ’ ) 

o 


where 


[«  - Sl),|  * c,  - if >V] 


VjW,  si  ♦.*  [«  - , <n  - Jf),eJ] 


The  program  for  the  calculation  of  the  exact  form  of  the  scattered 


field  is  now  formally  complete,  and  proceeds  by  two  stages.  First,  the  surface 


currents  are  calculated  from  eq.  (15).  Then  the  scattered  field  components  are 


calculated  at  an  arbitrary  field  point  by  eqs.  (20)  and  (23).  Both  steps  present 


formidable  computational  difficulties,  and,  for  the  case  of  an  arbitrary  surface 


T(£,n)  it  is  necessary  to  resort  to  numerical  analysis  for  usable  results. 


Before  turning  to  the  necessary  numerical  procedures  it  is  of  some  interest  to 


see  what  information  can  be  extracted  from  the  basic  theory. 


5.  SOLUTION  FOR  THE  CASE  [e  <<  l] 


As  mentioned  before  simple  iteration  of  eq.  (15)  does  not  produce  a 


power  series  expansion  of  iMC/h)  in  e.  However,  such  a series  expansion  can  be 


obtained  from  eq.  (15)  with  some  restrictions  on  T(£,n).  The  incident  field  is 


easily  expanded  in  powers  of  e, 


<t>  (C)  = eia^  COS  6°  [1  - iaET(C)  sin  eQ  - ya2£2T2(£)  sin2  dQ  ...] 


It  is  assumed  that  41  = i|/ ^ + ei|/ ^ + e2^2^  + ...  The  difficulty  arises  in 

the  expansion  of  H.!^  (ap  )/(ap  )3//2.  If  the  surface  satisfies  the  Lipshitz 
3/2  o o 

condition, 


It  - t • I < [(£  - £’  )2  + (n  - n')2]15 


then  (ap  )/(ap  )3//2  can  be  expanded  using  the  Gegenbauer  addition  formula, 

3/2  o o 


H3/2  (CtPo) 


^f-372  + 

(ap 


2 2 

where  p1  =[(£-£;')  + (n  - n')  ] . If  we  restrict  our  efforts  to  a second  order 

calculation  then  the  presence  of  a factor  of  e in  eq.  (15)  before  the  integral 
renders  the  0(e2)  term  effectively  proportional  to  e3  and  is  hence  discardable. 
Moreover,  M(C,#n';  £,n)  may  be  replaced  by  Mo(?',n'i  £,n)  for  the  same  reason. 

The  resulting  integral  equation  is  considerably  more  tractable.  Substitution  of 
the  assumed  power  series  expansions  in  £ yield  a set  of  equations  giving 

in  terms  of  i);  ^ , etc.  In  order  to  carry  out  the  integrations,  it  is 
necessary  to  introduce  the  Fourier  expansion  of  the  surface. 


T(£,n)  * / d<  / dp  c(K,y)  ei^K^  + Un* 


The  integrations  encountered  are  somewhat  complicated  and  much  use  is 
made  of  the  representation. 


ia  (PS+Xh) 


UPPPQ 


In  both  eq.  (29)  and  eq.  (30),  1 stands  for  the  unit  2x2  matrix.  The  zeroth 

order  term  ij/0^  - 26^"°*’  COS  9o  ^ *Xt  is  recognized  as  the  flat  surface  current  and 

o 

substitution  of  this  term  into  eqs.  (20)  and  (23)  gives  a reflected  wave  exhibiting 
the  oldest  known  fact  about  reflection,  angle  of  incidence  equals  angle  of 
reflection. 

From  eqs.  (20)  and  (23)  it  is  seen  that  the  e expansion  of  the  fields 

involves  the  e expansion  of  H (ctp  ')/(ap  ,)2/'2.  This  can  be  made  with  the  aid 

3/2  o O 

of  the  Gegenbauer  addition  formula  again  as  follows.  Introducing, 


pi 


[(5_^,2+  (n-^>2+  (if-)2]** 


(31) 


eq.  (21)  can  be  written 


p'.[Pl2  - 2p-  . (c,)2]'1 


(32) 


lez  * 

The  requirement  > et  (i.e.,  z'  > f ) allows  an  expansion  of 

^ a max  max 

(ap' )/(ap' )3/^2  in  powers  of  e,  again  using  the  Gegenbauer  addition  theorem. 

3/  2 o o 

This  expansion  is  carried  to  second  order  and  the  result  inserted  in  eqs.  (20) 

and  (23).  Moreover,  as  eq.  (22)  shows,  M®  (x’;  £,n)  is  a power  series  in  e to 

second  order,  as  are  Vx  and  V2  (eq.  (24)).  Using  the  power  series  expansion  in 

e for  4>(£,n)  results  in  a lengthy  calculation  which  determines  the  scattered  field 

as  a power  series  in  e.  Since  the  expressions  which  result  are  somewhat  cumbersome 

they  are  included  as  Appendix  I.  The  points  to  note  concerning  this  expression 

for  the  scattered  fields  are  (1)  the  result  is  valid  for  the  Lipshitz  condition 

of  eq.  (26),  and  (2)  the  result  is  valid  only  in  the  region  z'  > f x*  Of  course, 

2 

it  would  be  expected  that  the  solution  failed  if  e were  of  significant  size. 
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Several  other  features  of  these  e power  series  solutions  are  noteworthy. 


First,  with  the  assumption  of  periodicity  in  the  surface  so  that 


t (£,n) 


00  00 

l l 

m=_<»  n=-oc. 


iy-(m£  + nn) 
c e ii 
ran 


or,  equivalently. 


c(k,u) 


l 


n=-°° 


c 

mn 


6 (k 


(33) 


(34) 


the  expansions  in  e simply  reproduce  the  physical  optics  theory  of  Rice.  It  is, 

of  course,  not  clear  that  this  equivalence  persists  to  higher  order  than  two  since 

for  cubic,  or  higher,  terms  the  power  series  expansions  involve  corrections  to  the 

Green's  function  itself.  A second  result  is  implicit  in  the  £ expansion  of  the 
(1)  3/2 

Green's  function  H,  (ap')/(ap')  . The  result  is  a series  of  Green's  functions 

3/2  o o 

dependent  on  the  argument  p^.  It  can  be  seen  from  the  structure  of  p^  (eq.  (31)) 

that  the  series  is  interpretable  as  several  partial  Green's  functions  assembling 

the  field  from  a current  distribution  given  on  z'  = 0 = f . . This  is  true  in 

min 

spite  of  the  restriction  z'  > for  the  validity  of  the  e expansion  of  the  fields. 

When  the  condition  0 < z'  < f applies,  as  for  example  when  the  field 

max 

point  is  located  in  a "valley"  of  the  surface,  then  the  Gegenbauer  expansion  of 
H3/2  (^/(“Pq)372  requires  modification.  Consider  p^  as  a function  of  £ and  n with 
fixed  x',  y',  z'.  A cross-section  of  this  function,  together  with  the  surface  is 
shown  in  Fig.  2.  The  minimum  of  pj(£,n)  occurs  at  £ = kx'/a,  n = ky'/a  and  is 
equal  to  kz'/a.  For  large  £ and  n,  p|(£,n)  is  asymptotic  to  a cone  of  semi-vertex 
angle  ir/4  with  axis  £ = kx'/a,  n = ky'/a  and  vertex  located  at  z'  * 0.  In  general 
the  surface  defined  by  p^(£,n)  intersects  x(£,n)  in  a finite  number  of  regions 


* 

i '*♦ 


[ ■ 


near  £ = kx'/a,  h * ky'/a.  For  those  portions  of  the  surface  where  pj  < ex 
there  corresponds  a set,  denoted  by  y,  of  disjoint  regions  in  the  £,n  plane. 
For  the  regions  in  the  £,n  plane  outside  y,  which  is  denoted  y,  the  inequality 

n\  3/2 

p'  > ex  is  still  valid.  In  making  the  e expansion  of  H ' ‘ (ap’)/(ap')  for 
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region  y the  Gegenbauer  formula  is  used  exactly  as  before.  For  the  region  a 
new  length  is  introduced: 


p2  = C (£  ~ ~)2  + (n  - *£-)2  + (ex - 


max 


(35) 


so  that  p^  may  be  written, 


P’  - C(P')2  + 2p ' e (x  - x)  (—  - ex  ) + e2(x  - x)2]^ 

O z / max  a * n ma  v ma  v 


P 2 “ 


max 


max 


(36) 


Since  e(x  - x)  ^ ex  - kz'/a  the  Gegenbauer  addition  formula  may  now  be  used 
max  max 

to  generate  a series  of  partial  Green's  functions  whose  argument  is  ap^.  From  the 
structure  of  p^  it  is  seen  that  the  partial  Green's  functions  are,  in  the  region 
y propagating  waves  ("assembling  the  field")  from  current  distributions  located 
on  z'  = f (i.e.,  kz'/a  = ex  ).  These  currents  are  above  the  field  point 

and  as  seen  at  field  point  constitute  downward  propagating  waves.  The  rather 
mathematical  requirements  of  the  Gegenbauer  addition  formula  are  thus  seen  to 
produce  the  missing  downward  propagating  waves. 


6.  RESULTS  FOR  A PERIODIC  SURFACE 

In  this  section  an  analysis  of  the  case  where  x(£,n)  = x(C)  = x(£  + 2) 
is  undertaken  to  provide  further  comparison  with  the  physical  optics  solution. 
The  imposition  of  one-dimensional  variation  of  the  surface  is  made  solely  in  the 

15. 





interest  of  shortening  the  analysis.  The  techniques  to  be  presented  work  as  well 
for  the  full  case  t « x(£,n)  as  the  "one -dimensional"  surface  T ■ x ( 4 ) . The  choice 
of  period  equal  to  two  is  made  largely  in  the  interest  of  facilitating  the 
numerical  solution  of  the  source  integral  equation  (eq.  (15)).  Again  the  techniques 
presented  are  not  dependent  on  this  choice  of  periodicity;  only  numerical  quadrature 
formulae  need  be  altered. 

The  first  task  is  reduction  of  eq.  (15)  in  consequence  of  x = x(£). 
M(C',n';  £,n)  is  found  to  simplify  enough  to  eliminate  in  the  equation  for 
so  that  only  (|>  is  involved.  Since  ^(C,  n+A)  - ^2(C,n)  satisfies  a homogeneous 
integral  equation  (Fredholm  integral  equation  of  the  first  kind)  for  arbitrary  A 
it  can  be  assumed  that  <J>2(C,fl)  = ^(C).  This  proves  to  be  a sufficient  condition 
to  eliminate  \|»  in  the  equation  for  The  two  equations  thus  become  fully  de- 

coupled when  x = x (£) . All  that  is  necessary  for  reduction  to  two  Fredholm 
equations  of  the  second  kind  for  the  two  functions  (C)  and  >J>2(£)  is  the  71  integral 

of  Hi1,!,  (up  ) / (otp  )3^2.  This  is  easily  accomplished  using  the  contour  integral 
3/2  O O 

representation  of  eq.  (28).  There  results, 


^1(5’)  = 2ae 


ia[C’  cos  0 - ex'  sin  0 


1 ™ r H(1)(ap,)  , , 

J - i-j-  jd£  - - — (x-x')(l+e  xp  1^(5) 


*2(6')  =*  2beia^‘  cos  9o  " eT'  sin  9o-*  - i^-  jdf-  --1-  ^ --  [x-x'  - U-5')xc]  (4) 

* 1 (37) 

a and  b again  depend  on  incident  polarization,  and  with  x'  * x(£’}. 


P1  - [(£  - £')2  + e2(x  - x*)2]15 


(38) 


For  the  special  case  of  x(£)  « cos  «£  these  equations  have  been  studied  by  Zaki 


and  Neureuther. 


Either  of  the  eqs.  (37)  can  be  used  to  illustrate  the  method  of 
solution.  The  TM  mode  is  selected  here  and  accordingly  b * -1  so  the  relevant 
equation  is 


<|>U’)  = -2e 


ia[S ' cos  0 -ex'  sin  0 


H^W) 


-i  c f n iap.j 

J-  i~  d£  - - — [x-t'  - (C-C')te]*(C) 


where  the  2 subscript  has  been  dropped.  For  this  case,  since  a = 0,  4^(0  = 0. 
Since  4>(5  + 2)  = ip  (£)  it  follows  that  4*  (£  + 2n)  = 4 (Z)  for  n = 0,  ±1,  ±2,  ...  . 
Eq.  (39)  is  simply  a Fredholm  equation  of  the  second  kind  with  kernel, 

0 H.U)(ap.) 

k(C\5)  - -i  Y -■  - [t  - t*  - (C  - £')  tc]  (40) 


Because  of  the  periodicity  imposed  on  t(£)  the  kernel  is  seen  to  have  the 
periodicity. 


K (£ ' + 2n,  £ + 2n)  = k(£' » O 


As  it  stands  the  inhomogeneous  term  (in  the  sense  of  a Fredholm  equation)  has  no 

, , ..  . , . „ ia[£  cos  6^  - ct  sin  0 ] . . 

simple  symmetry.  However,  if  f(£)  * -2e  o o , then, 


f(C  + 2)  - el2a  COS  6o  f (?) 


so  the  restriction  of  0 to  satisfy  2a  cos  0 = 2irp,  p = 0,  1,  2,  ...  with 

o o 

irp  < a yields  the  same  basic  symmetry  as  the  surface.  Eq.  (39)  is  now  written 
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«K£*)  * f (C * ) + d£  K(£',£)  *(£) 


Fran  which  it  is  easy  to  show  that  <M£  + 2n)  - iM£)  satisfies  a homogeneous  integral 
equation  and  is  thus  equal  to  zero.  If  £'  is  confined  to  the  range  -1  < £'  < 1 and 
the  integral  written 


d£  K(£  ' , £ ) *{£)  ml  d£  K(£ ' , £ ) i|*  (£) 


£ + 2n)  ] i M£) 


*2.  00 

d£  [[  K(£', 
^ , n=-“ 


then  eq.  (43)  is  reduced  to  a form  that  is  suitable  for  numerical  solution. 


iM£')  = f (£ ' ) + d£  Kf(£\£)  «!»(£) 


where  the  folded  kernel,  Kf  (£',£)  has  been  introduced. 


K,(£',£)  = [ K(£',  £ + 2n) 


From  this  form  it  follows  at  once  that  Kf(£'  + 2n,  £)  = Kf(£',  £ = 2n)  = Kf  (£',£). 
It  is,  of  course,  possible  to  determine  Kf (£',£)  directly  from  eq.  (44)  but  the 
structure  of  K(£',  £ + 2n)  makes  the  convergence  of  the  series  very  slow  indeed. 

The  reason  for  this  may  be  seen  from  eq.  (40) . As  increases  so  does  the  term 
in  square  brackets  in  eq.  (40).  Therefore,  the  decrease  in  successive  terms  depends 
on  the  decrease  in  (apj)  which  is  proportional  to  p^,  a slowly  decreasing 

function. 


r ’ 


J 


It  proves  possible  to  improve  the  situation  by  using  the  integral 
representation,  eq.  (28),  together  with  the  Poisson  Sum  Formula.^  From, 


K (£',£)  - K(£',£)  + l [k(£\  £ + 2n)  + K(£\  £ - 2n)  ] (47) 

n=l 


and  eq.  (28)  there  results, 


l [K(£\  £ + 2n)  + K(£',  £ - 2n) ] = 
n=l 


2 2 2 
a 2 “>  n a 


J “ 2 [(£-£')tc  - (x-T ■ ) ] ~7-  J dv  V 2eV  4v  P1  £ e v cosh  (£-£') 


lcct 


n=l 


na 

v 


a 2 « 

I 

n=l 


2 2 
n a 


-2iea2T.  -r^r  / dv  v 2eV  4v  P1  J n e v sinh  (£-£') 

£ 2ni 


na 

v 


(48) 


Applying  the  Poisson  Sum  Formula  to  the  first  sum  on  the  right  of  this  equation  gives: 


2 2 2 _ 2 2 
00  n a 2 /—  , a . 2 00  , it  n 

V . na  , . /rr  ^ 7—  (£-£')  r i7rn(£-£  ) 5—  v 

lev  cosh  — (£-£')  = ->5  + — v e4v  Z e a2 

n=l  n=-°° 


(49) 


* 
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The  right  hand  side  of  this  substituted  in  the  first  integral  on  the  right  of 

2 2 

eq.  (48)  coupled  with  a change  of  variable  in  the  integral  v(l  - irn  /a  ) ■+•  v 
gives,  after  vise  of  some  basic  definitions  of  the  half-integral  Hankel  function, 
the  contribution, 

-K(C,?>  + *i[(E-{')t  - i — ! .lnl<5-«')+lat|T-T,|(l-  2T_)I’ 


t-t'  n=-°° 


19. 


(50) 


A 

A 


i"n(£-£')  + iaed-nV/cx^lT-T 


improvement  in  calculational  time.  It  should  be  noted  that  (£',£')  is  not  defined 
by  the  above  equation,  but  a simple  argument  based  on  the  fundamental  definition 


£ ) gives  K {£ 


7.  CONCLUSION 


An  exact  solution  for  scattering  from  a rough  surface  has  been  presented 


which  has  as  its  single  restriction  that  the  surface  is  perfectly  conducting. 
Calculation  of  the  surface  currents  forms  the  first  step  while  calculation  of  the 
scattered  field  from  these  as  sources  constitutes  the  second  step.  The  problem 
of  downward  propagating  waves  near  the  surface  has  been  shown  to  be  a natural 
consequence  of  step  two.  The  general  problem  was  then  limited  to  a surface 


exhibiting  periodicity,  f(x,y)  ■ f(x)  ■ f(x  + 2A)  and  improved  methods  for 
calculating  the  appropriate  kernels  of  the  source  integral  equations  were  pre- 
sented. The  results  of  the  theory  presented  were  shown  to  be  in  accord  with  the 
small  roughness  (e  <<  1)  theory  of  Rice. 
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APPENDIX  I 

Equations  (20)  through  (24)  give  the  scattered  magnetic  field  in  terms 

g 

of  surface  current.  From  eq.  (22)  the  e power  series  expansion  of  M (x',£,n)  is 

obtained,  as  is  the  e power  series  expansion  of  V^(x' ,Z,n)  and  V2(x’,C,n)  from 

eq.  (24).  The  e power  series  for  <M£,n)  has  been  given  in  eq.  (29).  The  remaining 

problem  is  the  e series  for  (ap ' )/ (ap ' ) 3/^2  with  p'  as  in  eq.  (21).  For  com- 

3/2  o o o 

parison  with  the  results  of  Rice  the  assumption  z'  > r is  made  and  the  Gegenbauer 

max 

addition  theorem  gives. 


"y; 

(ctp  • )3/2 
O 


23/2  r(3/2)  l (3/2  + £) 
1=0 


**1+3/2  (api}  JM-3/2(eaT)  (3/2)  kz1 . 

, . , 3/2  , ,£+3/2  C£  ap ' ' 

(ap|)  (eat)  1 


where  p^  is  given  by  eq.  (31) , and  c (x)  are  the  usual  ultraspherical 
polynomials:  c^2\x)  = 1,  c^2^  (x)  = 3x,  c^2^  (x)  = 3/2(5x2-l)  ...  etc.  3he 


e expansion  is  made  by  expanding  J^+3/2 


(eaT)/(eax) 


£+3/2 


and  retaining  terms  to 


order  e . We  set, 


4"  “v 

(apt>) 3/2 

o 


K + £K,  + C X, 
o 1 2 


Using  the  expansions  for  iMS,n)  and  MS(x’,£,n)  the  integrand  of  eq.  (20)  or  (23) 
can  be  written  as  a power  series  in  e where  each  integral  is  done  using  the  result 


derived  from  eq.  (28) 


- - “'“[(eV  V)11] 

*»*>•«*.>  ’ / « / TXTT^r  ‘1<P5nn) 

—oo  -oo  (£  +n  + p ) 


I 


(1-3) 


- 2*/f~l  (1  - p2  - XV/2'1/2  [6(1  - P2  - X2)1*] 

To  order  e2  ten  integrals  are  involved  for  (x1)  while  seven  suffice  for  (x') 
As  the  calculations  are  elementary  though  very  tedious  only  the  results  are  given. 


E(x,  y,  z;  p,  x)  = e 


ik[xp  + yx  + z(l  - p2  - X2)1*] 


s . . „ . ext  ik (x  cos  0O  + z sin  5o* 

% (X)  = Vo  e 


00  00 


(x)  = -ic2HoaJ  / dp  / dx  c (a  (p  - cos  8q) , a\)  E(x,  y,  z;  p,  x> 
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, 2 2 2 
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H®  (x)  = -e22H  a6  / dp  / dx  / dp'/  dx ' c( a(p-p'),  a(x~x')) 
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l-p  -x  + x(x  - x') 
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-p(x  - x') 
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